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Abstract. Observations of molecular gas in galaxies at all redshifts are critical for measuring the cosmic evolution in 
molecular gas density (Sl_f/2(z)) and understanding the star-formation history of the Universe. The ^^CO molecule 
{vrest for the J=l— >-0 rotational transition = 115.27 GHz) is the best proxy for extragalactic molecular hydrogen 
(H2), which is the gas reservoir from which star formation occurs. The detection of CO at high-redshift has been a 
' growing industry over the last few years with successful molecular gas detections out to z^Q. Typically, redshifted 

» I , high- J transitions are observed using mm-wave instruments, with the most commonly targeted systems being those 

■ with high star formation rates such as submm galaxies, and far-infrared-bright quasars. While the most luminous 
objects are naturally the most readily observed, observations of the typical members of the galaxy population 
which exhibit modest star-formation rates are essential for completing the picture. The arrival of ALMA will be 
revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include more normal 
star-forming galaxies, however the limited field-of-view when observing at such high frequencies makes it difficult 
to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces 
a strategy for a systematic search for molecular gas during the Epoch of Reionization (z'^7 and above), capitalizing 
on the fact that the J=1^0 transition of ^^CO enters the upper frequency bands of cm-wave instruments at the 
appropriate redshift. The field-of-view advantage gained by observing at such frequencies, coupled with modern 
broadband correlators allows significant cosmological volumes to be surveyed on reasonable timescales. In this article 
we present an overview of our proposed observing programme which has been awarded 6,500 hours as one of the 
Large Survey Projects for MeerKAT, the forthcoming South African Square Kilometre Array pathfinder instrument. 
Its large field of view and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such 
work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility 
of such observations and predict detection rates. 



1. Introduction 

Carbon monoxide (CO) is the most abundant molecule 
after molecular hydrogen (H2), however the prohibitively 
high excitation temperature of the latter means that it 
is generally unobservable except in hot, shocked environ- 
ments (e.g. Rigopoulou, 2002). As such the CO molecule 
is generally used as a proxy tracer for H2, with a CO/H2 
ratio within the Milky Way of '^10"''. Observations of 
the CO molecule are thus critical in order to determine 
the cosmic evolution of the molecular gas density of the 
Universe (r^jjj (z)). Since the molecular gas in a galaxy is 
the fuel reservoir for star-formation such observations are 
vital tools for studies of galaxy formation and evolution, 
and for tracing the star-formation history of the Universe. 
See the article by Walter et al. (2011) for a recent review 
of the importance of detecting molecular gas at high-z. 

The usual targets for such observations are the 
rotational-vibrational transitions of the ^^CO molecule. 
The J=1^0 ground-state transition has a rest-frequency 
of 115.27 GHz and higher-J transitions occupy integer 
multiples of this frequency. As such, high-redshift obser- 
vations are generally performed with millimetre-wave in- 
terferometers, targeting J — >• (J-1) for J>1 (e.g. Tacconi 
et al., 2008), although the advent of the Expanded Very 
Large Array (EVLA) with its new Ka-band receivers (26 - 
40 GHz) has resulted in numerous recent detections of 
the ground-state transition of ^^CO at high-z (e.g. at 
z~2, Ivison et al., 2010; z^3, Riechers et al., 2010). 
Observations of the ground-state transition provide the 
best estimate of the total amount of molecular gas within 
a galaxy. Typical targets for high-z observations to date 
have either been objects which are undergoing episodes of 
vigorous star-formation such as "sub-mm (selected) galax- 
ies" (e.g. Blain et al., 2002) with star- formation rates in 
excess of 1000 Mq yr~^, or luminous quasar host galaxies 
(e.g. Riechers et al., 2008). 

It is desirable to observe multiple CO transitions in 
a given object in order to accurately determine the gas 
excitation state. Coupling this with the need to perform 
such observations over a range of redshifts across cosmic 
time implies an observing frequency range which essen- 
tially spans two orders of magnitude, from local, high-J 
observations at close to 1 THz, to ground-state observa- 
tions at 10 GHz which target objects at Z'^10, well within 
the cosmic Epoch of Reionization (EoR; e.g. Dunkley et 
al., 2009). 

The MESMEbQ survey is one of the approved Large 
Survey projects for the forthcoming MccrKAT radio tele- 
scope, granted 6,500 hours of observing time, with the 
primary goal of detecting large numbers of C0(J=1— t-O) 
emitters during the EoR. This article presents an overview 
of this observing programme, and couples an existing sim- 
ulation of the evolving molecular line properties of high- 
redshift galaxies with some instrumental considerations in 
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order to demonstrate the feasibility and predict the return 
of such observations. 

2. MeerKAT 

The MccrKAlU array is currently under development 
and construction in the Karoo region of South Africa's 
Northern Cape province. It is one of the two complemen- 
tary dish-based pathfinder instruments under construc- 
tion by the two potential host nations of the core of 
the Square Kilometre Array (SKA), the other being the 
Australian SKA Pathfinder (ASKAP). 

The current specification of MeerKAT describes 64 
parabolic antennas with a diameter of 13.5 metres, and 
offset Gregorian optics. The unblocked aperture offered 
by such a system results in improved primary beam per- 
formance, and this optical configuration is a recent shift 
in design from having receivers at the prime focus, as 
is the case for KAT-7, the currently operational seven- 
element MeerKAT prototype. When complete, MeerKAT 
will be the most sensitive centimetre-wave instrument in 
the Southern Hemisphere, comparable in sensitivity to the 
EVLA, and much faster in mapping speed by virtue of its 
smaller dishes and hence larger field-of-view. 

Final observing frequency specifications are likely to 
consist of three receiver bands: 0.58 - 1.015 GHz, 1 - 1.75 
GHz and 8 - 14.5 GHz, with a correlator capable of de- 
livering complete coverage of the two lower bands and up 
to 4 GHz of instantaneous bandwidth for the high-band, 
split over up to 25,200 frequency channels. 

The proposed dish layout results in a fairly compact 
core, with over 50% of the dishes within a 500 metre diam- 
eter and 75% within a kilometre, and a maximum baseline 
of approximately 8 kilometres. A spur of antennas run- 
ning along the access road (in an approximately easterly 
direction from the core) is also under consideration which 
would boost the maximum baseline to ~50 kilometres. 

3. The S^-SAX simulation 

The S^-SAX simulation (Obreschkow et al., 2009a and 
references therein) is built upon the evolving dark mat- 
ter skeleton of the Millennium Simulation (Springel et al., 
2005), and the galaxy semi-analytics of De Lucia & Blaizot 
(2007), whereby the evolution of model galaxies placed at 
the centre of the dark matter haloes was tracked. A physi- 
cal recipe was applied to these two existing layers in order 
to divide the cold gas in the model galaxies into H, H2 
and He, and realistic gas distributions and velocity pro- 
files were also added to the H and H2. 

Obreschkow et al. (2009b) subsequently applied a 
model to predict the luminosities of the first ten J transi- 
tions of CO, including molecular cloud geometry and over- 
lap, and contributions to gas heating from active galactic 
nuclei and starburst activity, as well as the cosmic mi- 
crowave background which becomes increasingly impor- 
tant at higher redshifts. 

^ http : / / www ■ ska . ac . za| 
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Fig. 1. The velocity-integrated C0(J=1— t-O) line flux for ~10,000 simulated galaxies as a function of redshift. Key 
predictions of the S'^-SAX simulation are the evolution with redshift of the line strength and the ratio of atomic to 
molecular gas (the latter represented by the colour gradient), and the cosmological decline in the source counts. The 
redshift ranges of some key spectral lines are plotted, as probed by an observation which covers 10.5 - 14.5 GHz. Note 
that the vertical positions of these markers has no significance. 



The simulation results in ^30 million galaxies at z = 0, 
each of which has a well-defined merger history trace- 
able to very high redshift, and a mass limit approximately 
equal to that of the Small Magellanic Cloud. The database 
is publicly- accessiblelfl and contains amongst other things 
the cosmologically evolving atomic and molecular line 
properties of every simulated galaxy. 

Analytic expressions can be used to derive model line 
profiles from parameters within the database, and adding 
appropriate instrumental noise to the line profile can allow 
us to make some predictions about the detectability of the 
line. This approach is employed in order to obtain a first- 
order prediction for the return of the MESMER survey. 
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4. MESMER 

The primary objective of the MESMER survey is to em- 
ploy the high-frequency (14.5 GHz and below) receivers of 
MeerKAT to perform a systematic search for emitters of 
the J=1^0 line of ^^CO {vrest = 115.27 GHz). 

The upper limit of the high-frequency band at 14.5 
GHz corresponds to C0(J=1— >0) at z = 6.95 which is 
inside the tail-end of the Epoch of Reionization. Using the 
wide bandwidth of MeerKAT to extend the observations 
to lower frequencies pushes the targets further into the 
EoR. An observation with 4 GHz of bandwidth between 
10.5 and 14.5 GHz would target CO{J=1^0) in the range 
6.95 < z < 9.98, as delineated on Figure [1] 

The efficiency of a given survey is defined by the vol- 
ume of space that can be covered and the fiux limit that 
can be reached during the allotted time. Figure [5] demon- 
strates the approximate size of the primary beam, full- 
width at half-maximum, for three interferometers target- 
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Fig. 2. Approximate fields-of-view (full- width at half- 
maximum) for the ALMA, EVLA and MeerKAT inter- 
ferometers, tuned to frequencies which target the CO 
molecule at zr^l. The advantage of using cm-wave in- 
struments for efficient line searches at this redshift should 
be immediately apparent. Also plotted within the squares 
for reference are the Hubble Ultra-Deep Field (Beckwith, 
2006) and the Hubble Space Telescope observations of 
AbeU 1689 (Benitez, 2002) 



ing CO at z^7. Observing at the lowest feasible frequency 
maximises the field of view, and to target CO at Z'-^T, 
ALMA would be sensitive to the J=7— )-6 transition at 
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~100 GHz. Figure [2] demonstrates the advantages that 
MeerKAT has for covering large sky areas for observations 
of this nature. Also shown for reference are the Hubble 
Ultra-Deep Field (Beckwith, 2006) and the Hubble Space 
Telescope observations of Abell 1689 (Benitez, 2002). 

As is common practice for many survey programmes, 
a tiered or "wedding cake" approach is adopted, whereby 
there are a few distinct observing strategies, and a trade- 
off is made between sky coverage and survey depth be- 
tween them. The three tiers are described in Section [S] 

5. Observing tiers and predicted detections 

This section describes the current plan for the observing 
strategy, and how the specifications of MeerKAT can be 
combined with the S^-SAX simulation in order to predict 
the yield of a given observation. 

Each pointing (or set of pointings) probes a volume 
of space defined by the field-of-view of the instrument (or 
size of the mosaic) in the spatial directions, and by the 
bandwidth of the observation in the redshift direction. A 
sub-catalogue of galaxies within an equivalent volume of 
space defined by a particular observation is then extracted 
from the S^-SAX database by submitting a SQL query via 
the web interface. 

The root-mean-square (RMS) noise in an image pro- 
duced from a single channel, with naturally-weighted 
visibilities, can be calculated via the standard sensitiv- 
ity equation (e.g. Equation 6.62, Thompson, Moran & 
Swenson, 2001): 



where fc^ is the Boltzmann constant, Tgys is the system 
temperature (50 K), A is the effective area of the 13.5- 
metre diameter antennas, tjq is an efficiency term due 
to correlator quantization losses (0.95), Up is the num- 
ber of averaged orthogonal polarization products (2), N 
is the number of antennas (64) Av is the channel width 
(3.8 MHz, i.e. 1024 channels for 4 GHz of total bandwidth) 
and T is the on-source time. 

The analytic line profiles for each galaxy in the sub- 
catalogue are then evaluated sequentially. A source is con- 
sidered to be detected if (a) the fine ffux integrated over 
the channel width at the peak of the line exceeds 5cr; and 
(b) the full- width of the line at the half- maximum level 
exceeds three channel widths. The threshold value of 
suggests that the resulting source counts will be conser- 
vative, however this may be offset by the assumption that 
the data can be perfectly calibrated. 

5.1. Tier 3: a square-degree survey 

The broadest and shallowest observing tier is likely to be 
a square-degree survey over a field with excellent multi- 
wavelength coverage (e.g. the Extended Chandra Deep- 
Field South). It is anticipated that 2,000 hours will be 
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Fig. 3. The velocity integrated C0(,/=1— >0) line fluxes 
plotted as a function of redshift from the simulated galax- 
ies within a volume of space defined by one square degree 
and 6.95 < z < 9.98. Detections are determined by the 
criteria outlined in Section [5] and these are marked in red. 
A histogram (with a logarithmic y-axis) of the sources in 
50 redshift bins is plotted at the top of the figure. 



dedicated to this. With a traditional hexagonal mosaic 
pattern and an assumed calibration overhead of 20%, this 
results in approximately 340 pointings with 4.7 hours of 
actual on-source time each. 

The channel RMS for such an observation is 32 /iJy, 
and the extracted square-degree of the S'^-SAX simulation 
at 6.95 < z < 9.98 (4 GHz bandwidth) contains 87,597 
galaxies. Using the method outlined above suggests that 
such an observation would detect 186 C0( J=l— )>0) emit- 
ters. Note that no primary beam effects are taken into 
account as the mosaic pattern should allow coverage of 
the field with approximately uniform sensitivity. 

Figure [3] shows the velocity-integrated C0( J=1^>0) 
line fluxes for the simulated sources as a function of red- 
shift. Detections are marked in red, and binning these de- 
tections by redshift forms the green histogram at the top 
of the figure. 

5.2. Tier 2: quasars and Lyman-a emitters 

Tier 2 will consist of 9 pointings, each having a total ob- 
serving time of 300 hours, which with the calibration over- 
head results in 240 hours of on-source time. The channel 
RMS in this case is 5.3 /iJy. Targets are likely to be the 
highest-redshift quasars, regions containing strong over- 
densities of Lyman-a emission, and regions selected for 
future, deep infrared observations with the James Webb 
Space Telescope. Although the current highest-redshift 
quasars are outside of the MeerKAT frequency range, 
Willott et al. (2010) using current constraints on the 
quasar luminosity function at z>6 predict that the VISTA 
survey programmes will detect >20 quasars at z>6.5, and 



Heywood et al.: MESMER 



>25% of these are expected to fall within the redshift 
range observable by MeerKAT. 

The most luminous quasars are likely to correspond 
to the most molecule-rich galaxies in the early Universe 
(Walter et al., 2009) and narrowband search techniques 
which discover overdense regions of Lyman-a emission are 
likely to be revealing the sites of primordial cluster forma- 
tion (Ouchi, 2005). 

It is likely that such deep observations will not only 
reveal the molecular gas of the primary target, but will 
also contain many other CO-linc emitters. Selecting sev- 
eral random volumes of space corresponding to a single 
pointing with 4 GHz of bandwidth suggests that each 
pointing is likely to yield a total of ^20 such detections. 

The primary beam of an interferometer causes a drop 
in sensitivity away from the pointing centre and this is 
approximately accounted for when simulating this tier. A 
normalised Gaussian with a FWHM corresponding to the 
size of the Airy disk for a 13.5 m aperture is centred on the 
pointing direction. This then defines the function by which 
source line fluxes are attenuated based on their position 
in the field. 

5.3. Tier 1: strong gravitational lenses 

The deepest pointings will target three strong gravita- 
tional lenses. Potential targets are Abell 370, Abell 1689 
and MACS2129 which have no bright radio galaxies in 
the field. The magnification afforded by the foreground 
cluster will amplify the strength and counts of the back- 
ground high-redshift population. Each pointing will have 
600 hours of observing time (480 hours on-source) , giving 
a channel RMS of 3.8 /iJy. 

Mass models for the above galaxy clusters applied to 
the dN/dz counts extracted from the S^-SAX simulation 
suggest an average magnification factor of ~2 across the 
field, with high factors ( 5-30) in small areas of the source 
plane. It is estimated that the overall detection rate for 
high-z CO emitters for these observations will be ~30-40 
galaxies per pointing. 

In the small regions of very high magnification galaxies 
with integrated line fluxes of '^2xl0~^ Jy km corre- 
sponding to objects with H2 masses as low as ~10^ M0 
and very modest levels of star formation (~5 Mq yr^^) 
start to become detectable. The boosting of fluxes at the 
higher end of the luminosity function emphasises the im- 
portance of targeting several lens flelds to account for cos- 
mic variance. 

Note also (with reference to Figure [2]) that the fleld of 
view of MeerKAT is sufficient to cover the whole strongly- 
lensed area of such clusters in one shot. There is no need 
to divide the observing time over multiple pointings with 
painstaking targeting of the lens caustics. 

6. Parallel science goals 

Very deep, high-frequency observations as will be gener- 
ated by this programme will obviously have a wealth of 



other applications. In addition to the primary goal of de- 
tecting molecular line emission during the EoR there will 
also be the possibility to search for other spectral lines at 
lower redshifts. 

Hydrogen cyanide (HCN( J=1^0), i^rest 88.63 GHz) 
and formylium (HCO+( J==1^0), Vrest = 89.19 GHz) lines 
will be visible for the redshift range 5.1 < z < 7.5, cor- 
responding to 14.5-10.5 GHz. These lines are tracers of 
dense (nj/^ > 10^^-'^"^^) m~^) regions of gas, which are 
typically associated with the star-forming cores of giant 
molecular clouds (see e.g. Solomon and Vanden Bout, 
2005), and the ratio of these lines is often used as an in- 
dicator of starburst versus AGN activity in galaxies (e.g. 
Imanishi et al., 2007; Knudsen et al., 2007). 

Despite these lines being significantly weaker than CO 
by a factor of 10-40, predictions from Carilli & Blain 
(2002) show that a small fraction of the detected high- 
z C0( J=l— 7-0) lines will be false detections corresponding 
to HCN at lower-z. Those authors estimate a detection 
rate for HCN (expressed as a percentage of CO detections) 
of 1.5% at z = 1, rising to 5% at z = 3. 

Separating these line sources from the genuine z>7 CO 
detections can be accomplished by identifying both HCN 
and HCO^ in the same spectrum due to their proximity 
in frequency. For detections where only a single line is 
visible then near-infrared identifications and photometric 
redshifts are likely to be necessary. Failing that, follow-up 
observations (e.g. with ALMA) will confirm the redshift 
by targeting other spectral lines. 

Additionally the data will facilitate a deep search for 
0.5 < z < 1.2 water masers {h'rest — 22 GHz). Currently, 
the two most distant detections are at Z'^2.5 and z~0.66 
(Impellizzeri et al. 2008; Barvainis & Antonucci 2005) the 
former being gravitationally lensed by a factor /i ~ 33. 
The redshift range probed by MESMER provides an am- 
bitious yet practical next step to push the H2O luminos- 
ity function out to higher z. Any detections will be fol- 
lowed up with VLBI to measure robust kinematic black 
hole masses and independent geometric distances, pro- 
vided the maser emission is sufficiently luminous when re- 
solved. Independent distance measurement of course pro- 
vide an estimate of the Hubble parameter in a cosmologi- 
cally significant redshift range where dark energy becomes 
comparable to the matter energy-density in the Universe. 

The redshift ranges for the two parallel spectral line 
surveys described above are marked on Figure [TJ 

The MESMER programme will also yield very deep 
~12 GHz continuum data. Spectral indices determined by 
combining these data with those of lower-frequency ra- 
dio observations (e.g. from the MIGHTEE Large Survey 
project. Pis van der Heyden & Jarvis) will be very useful 
for separating the synchrotron from the free-free emission 
in star-forming galaxies, and for probing AGN / starburst 
composites. 
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7. Imaging simulation pipeline 

The next step in terms of simulating the observations de- 
scribed in this document is to move from the simple ana- 
lytic approach to a full imaging simulation, and a software 
pipeline has been assembled for this purpose. The pipeline 
takes as its input a subcatalogue of galaxies from S^- 
SAX and a list of observational parameters, detailing such 
things as the observation duration, the pointing direction 
and the frequency setup. The basic steps in the simula- 
tion chain are as follows: (a) generation of an "ideal" ra- 
dio sky from the simulation database using the S3Tools 
(Levrier et al., 2009); (b) pixelwise attenuation of the sky 
model by an analytic primary beam model; (c) use of the 
CAs41 sm tool to generate a per-channel Measurement Set, 
to serve as an empty "crate" for visibilities; (d) use of 
the MeqTrees package (Noordam & Smirnov, 2010) to in- 
vert the sky model and fill the Measurement Sets with 
appropriately- noisy visibilities; (e) imaging and (optional) 
deconvolution of the visibilities using the Iwimager, part 
of the casarest package; and (f) (optional) correction of 
the primary beam model. 

This process can operate from start to finish without 
user intervention using a combination of shell and Python 
scripts, and by operating MeqTrees in batch mode. Each 
frequency channel is treated as an individual simulation, 
which naturally lends itself to the jobs being split over 
multiple processors. The final set of channel maps can also 
be optionally combined into a single FITS cube. 

This software chain has been used to simulate a sin- 
gle 480-hour pointing with 2 GHz of bandwidth (12.5 - 
14.5 GHz), consistent with a MESMER tier 1 observation 
of high-z C0(J=1— >-0), with 2:1 channel averaging. Note 
however that the effects of the foreground lensing clusters 
are not (yet) included. The end result is a 4096 x 4096 
X 256 pixel datacube. The final FITS cube is 17 GB al- 
though the whole process generates approximately 0.3 TB 
of data, mostly visibilities. A solid detection from within 
this simulated cube, showing the C0(J=1— >-0) emission 
from a pair of galaxies at z = 7.02, is presented in Figure 

El 

The goal of this work, in addition to obtaining a more 
robust estimate of the detection rate, is to test the ef- 
ficacy of source finder software and spectral-line stacking 
code in the presence of known noise and sky properties for 
MeerKAT-scale datasets. A separate paper by Heywood et 
al. (2011; in prep) will describe this work in full. 

Note that although the simulation method is described 
here in the context of MeerKAT it can be used to sim- 
ulate completely arbitrary interferometers, and MeqTrees 
allows for the introduction of complex instrumental effects 
(e.g. more accurate primary beam treatment from electro- 
magnetic dish simulations) without having to develop code 
beyond the Python layer. Additionally, the sky model can 
also be completely arbitrary, but without modification to 
this collection of scripts, it can be generated from the HI 
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Fig. 4. A 180 X 180 x 7 subcube extracted from the full 
size single-pointing simulation, showing the detection of 
C0(J=1— i^O) emission from a pair of z = 7.02 simulated 
galaxies. Channel number and frequency are marked on 
each frame. The horizontal bar on each frame spans 5 
arcseconds. 



component of S^-SAX, or any of its other CO transitions, 
and also the S'^-SEX continuum simulatior(f| devised by 
Wilman et al. (2008). The article by Levrier et al. (2009) 
describes in more detail the possibilities for generating sky 
models from the S^ suite of simulations. 



* [http : //casa ■ nrao . edu] 
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8. Conclusions 

The MESMER Large Survey project will utilise the high- 
frequency receivers of the forthcoming MeerKAT array 
and, using a combination of very deep single pointings and 
shallower wide-area observations totalling 6,500 hours, is 
expected to detect the J=:1->G line of the ^^CO molecule 
in ^^400 galaxies in the range 6.95 < z < 9.98, during the 
cosmic Epoch of Reionization. 

The survey will naturally include parallel spectral line 
surveys at other redshifts, notably a HCN / HCO+ at 
5.1 < z < 7.5, and H2O masers at 0.5 < z < 1.2, as well 
as deep '^12 GHz continuum maps which will be valuable 
for AGN / starburst studies. 

Making predictions about the feasibility of such obser- 
vations (and thus getting a handle on the science return) in 
addition to investigating via simulation the instrumental 
performance of current and next-generation interferome- 
ters was the primary motivation for the development of 
the S'^-SAX simulation (and the other S"^ simulations). 

This article has presented both a simple, first-order ap- 
plication of the simulation via the generation of analytic 
line profiles, and also a framework for extending the in- 
strumental simulations into a full interferometric imaging 
simulation. The latter approach is particularly powerful 
for assessing current instrumental designs, testing software 
on large-scale datasets, and naturally provides a more ro- 
bust prediction as to the high-z molecular line detections 
that can be made with the MeerKAT array. 

The framework used for this investigation can also be 
used to simulate other key science, including extragalactic 
neutral hydrogen and continuum. 
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